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Abstract. We continue a previous paper to show that Mel'nikov's 
first order formula for part of the separatrix splitting of a pendulum 
under fast quasi periodic forcing holds, in special examples, as an 
^ , asymptotic formula in the forcing rapidity. 

5— i ' §0. Introduction. 

Recently there has been renewed interest in a problem treated in the paper [G3] . One of 
the several questions posed in [G3] was to find upper and lower bounds on the splitting 
between the stable and unstable manifolds of the invariant torus with rotation vector 
u; = lo in the Hamiltonian system ( ll Thirring modeV): 
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m; V 2 u -A+—A-A+ — +g 2 J (cos^-l) + [if(a, V ), (0.1) 
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where (I,<p) 6 E x T, (A , a ) T are canonically conjugated variables, 

UqCR 1 1 , and J, Jo > (respectively, "rotators' moments of inertia" and "pendulum's 
moment of inertia"), g > (g 2 is the "gravity"); w ,/x are parameters. Here J could 
be a scalar or a diagonal (/ — 1) X (/ — 1) matrix. And setting v = (n, u) S 1 , \u\ = 

M + |£| = |n|+EtiN: 



oo 

On 

q . f(a,(p)= 22 fa,nCos(E ■ a +ncp) , (0.2) 

43 ' 

O ■ where / w are fixed constants and Nq > is fixed. 1 

The paper [G3] did solve completely the upper bound question for / a trigonometric 
polynomial: in the sense that it derived generically optimal bounds for the Fourier trans- 
form of the 2-dimensional splitting function A*(a) (a vector valued function defined in 
[G3], p. 366). 

Subsequently it was rightly pointed out, [DGJS], that it would be interesting to just have 
examples in which some lower bound could be computed as given by the naive first (non 
trivial) order perturbation theory prediction (i.e. "Mel'nikov integral"). For instance 
one could study the difference in the free pendulum energy ho = 2j~I 2 + Jog 2 (cos <p — 1) 
evaluated in two "corresponding points" of the stable and unstable manifolds of the 
invariant torus. Such difference is a very special case of the splitting functions considered 
in [G3]. 

Of course studying the splitting only via the variations of a single function ("observ- 
able") is very reductive. The stable and unstable manifolds are 2-dimensional surfaces, 
if observed on a (4-dimensional) section transverse to the motion and at fixed total en- 
ergy, hence it is clear that studying the difference of the two manifolds at "corresponding 
points" on the chosen section requires the simultaneous analysis of 2 observables. Thus in 



It seems clear that the often sought "first order of perturbation theory" dominance as rj — * can only 
hold if No < oo and fu ,n is nonzero and "quite large" for "many" u_ 's. Hence Nq < oo is a reasonable 
assumption, together with the assumption that we introduce later on f v . 
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[G3] two independent observables are considered (more are not necessary: the dimension 
of the manifolds being 2). 

One also needs to have good control over the Mel'nikov integral, i.e. on the usually 
explicitly known first order (in e) expression of the splitting: not so easy in general. 
Control is possible in some rather special cases like the 3 degrees of freedom system 
(i.e. two rotators) with f(a,ip) = J2 V fa ( cos E 1 a)(cosif), where /,, = e~ K l- or 
more generally F'e~ K l-l < \f„\ < Fe~ K l-l for some F',F,k > 0, and lu_ = (^01,^02) a 
rotation vector with a golden mean rotation number woi/^02 = ^(V5 — 1). 

The problem has been studied, [DGJS], without taking into account [G3] but without 
any saving of work because a theory of the upper bound equivalent to [G3] has to be de- 
rived one way or other: and in fact they had to restrict considerations to the isochronous 
cases (J = +00): taking J < +00 does not change the first order analysis but changes 
completely the higher orders treatment. 

We show here that the analysis of [G3] yields an upper bound on the splitting optimal 
and sufficient to deduce the dominance of Mel'nikov's approximation for the part of the 
splitting defined by the variation of the function ho above defined, if the perturbation is: 

f(a,tp)= ^ f„ cos u ■ a cos <p, N d = rf 1 , (0.3) 

v, \v\<N 

with f E > Fe~ K l-l for some F, n > and all v_? This differs from the function 
considered in [DGJS] by a quantity 0(e~ r ' 1r ) which is much smaller than any expected 
splitting. The latter is expected to be given by the Mel'nikov integral, hence to be of 

order > 0(e~ o ^~^). 

Extending the sum over v to all y_ 's would require extending the results in [G3] to the 
case in which the perturbation is analytic, while one of the "philosophical assumptions" 
of [G3] was that one should understand first only polynomial cases. However the purpose 
of [DGJS] was to provide an example of the Mel'nikov integral dominance and the (0.3) 
is as good for the purpose. 

Generalizing [G3] to non polynomial analytic cases is easy: but to limit the length and 
the simplicity of the present comment on [G3] we do not provide the details. We can in- 
dicate that an extension to analytic cases has been performed in the very similar problem 
studied in [BCG] (providing again an example of the Mel'nikov integral dominance). 

The fact that the bounds in [G3] are uniform in the matrix J (which is summarized 
in [G3] by calling twistless the invariant torus of (0.1) with rotation vector u) will 
imply, below, that the perturbation / in (0.3) also generates a splitting measured by 
the variation of the function ho that is asymptotically exactly given by Mel'nikov's first 
order perturbation result for all J > Jo- The treatment given below being more general 
than what would be needed to just obtain the splitting of ho also solves (in part: see 
below) the problem posed in [RW]: a paper in which the basic strategy is adapted from 
[G3] but which seems to contain incorrect usage of [G3] leading to errors in intermediate 
steps and in the final results, see [GGM3]. 

The following section assumes full knowledge of [G3] and we consider it just as a (con- 
ceptually trivial) final comment to it, made interesting by the idea of [DGJS] of studying 
lower bounds for single observables: therefore we call it §9 (§10 below contains a few com- 
ments to relate the results to intervened papers), as the paper [G3] ends rather abruptly 
(because of exhaustion of the author, and of the problem) at §8. We use the notations of 
[G3] and we freely quote the formulae in the eight sections of [G3] . The following section 
has simply to be regarded as a new section of [G3]. The above is an introduction and 
formulae have been labeled by to avoid confusion with the labels used in [G3] . 



This note is made necessary by recent claims that [G3] would not be sufficient to get an upper bound 
which would yield the full result of [DGJS], see [St]. 
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The relation between the splitting of the observable ho and the 2 dimensional splitting 
A|(a) introduced in [G3] is, on the Poincare section considered in [G3] (i.e. at <p = it, 
and fixed energy) and if J = +00, simply - w • A f (a). In the anisochronous case 
(J < +00) it is slightly more involved. Unless explicitly stated, we shall mean that the 
"splitting" is the vector At(a) of [G3], p. 366. 

§9. Homoclinic splitting in presence of several fast rotators. 

An interesting consequence of the theory of §8 is an expression of the splitting in the 
I = 3 case and when the vector to has all components fast, i.e. lo = r\~^ uj_ and rj — > 0, 
while w is fixed and verifies the Diophantine condition (1.3). 

Examining the tree expansion for A|(a) one realizes that A|„ is given by a sum of 

tree values each of which is bounded by £ = e~' — '- 9 (f~ d ) (where d is prefixed and 
chosen equal to 772, as in §8, to treat the case I = 2) with a in general different from 
v_ but nonzero, see p. 379. 

Each tree will in general contain several bubbles (see §6, p. 376). Here it will be 
convenient not to lump together the contributions from the trees with the same free 
structure: rather it will be convenient to leave them separate: thus a tree value will be 
described as in §6 but its fruits will contain their seed inside (see §6, p. 377). We should 
call them seeded fruits, for obvious reasons; below we shall simply call them just "fruits" : 
they can be either ripe or dry, as defined in §6, p. 376. Therefore we redecompose the 
resummation trees into single trees by specifying the seed inside each fruit. This is a 
minor change and it is convenient as it will show that the bound in (8.1) can be trivially 
and greatly improved in the case of perturbations like (0.3). 

The small factor £ arises for some y_' with \ v!\ < Nh (N being the cut-off parameter 
N = i] 1 ) m two different ways. 

(1) In trees with free momentum u' (sec §6, p. 375) and only ripe fruits (see §6, p. 376) 
it arises by bounding carefully (i.e. by complex integration) the integrals relative to the 
free nodes (called v) variables t v and by using the trivial bounds to bound the fruit values 
(i.e. with no excursion into the complex plane to estimate them, which of course would 
not be allowed because of lack of analyticity of the integrals that express the values of 
ripe fruits). 

(2) In trees with at least one dry fruit it arises from a fruit located inside it, perhaps 
"very deeply" (i.e. inside many other fruits), which only contains ripe fruits. It exists 
necessarily (note that no fruit at all is a special case of only ripe fruits) and in this 
case the free momentum cannot vanish (see §6): this was the reason why the induction 
described in §6 worked. 

Therefore there is a natural resummation of the splitting series ^ h e h e 1 -' — A j v : 
namely we separate the first order contribution and we collect together all contributions 
from trees whose value is bounded by the same small factor £. In this way: 

00 

A T (a)= E e^ e -£l^l( M A; £ +X> fe A^(a)) , (9.1) 

vjtO h=2 

where e'-'-e - ^'-'-'^^ isolates the first order terms, while the quantity e 1 -' — 

e"^ 1 -'- 1 A^(a) is defined as the sum of the values e^-'-e - ^'-'-' A ^ (a ) of all 
trees •& of order h containing cither only ripe fruits (or no fruits at all) or a dry fruit 
which has inside a seed that is a tree with only ripe fruits (or no fruits at all); see items 
(1) and (2) above. 

Since a tree can have several fruits which can be "exponentially bounded" there is am- 
biguity in attributing the terms which contain several small factors e" 5 ?'- '-j with mo- 
menta v_ x , v_ 2 i ■ ■ ■■ m the latter cases we just make an arbitrary choice to attribute such 
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terms to A^ (a) with some y_ among v_ 1 , v_ 2 , . . . (e.g. v is the first among u 1 , v_ 2l . . . 
in some lexicographic order). 

The series involving \i h A y v (a) starts at least at order two because the trees must be 
at least of second order, having already separated out the first order contribution. Fur- 
thermore each tree contains as a factor the product of the couplings fv that correspond 
to its nodes v with order label S v = 1 (see p. 367); and the set of the v_ v must be such 
that a suitable subset W of the set of nodes verifies J2 veW —v = This means that 
we can bound the contribution of each tree d with h nodes v of order 5 V = 1 (and hence 
total number of nodes m < 2h) by the bounds in Appendix Al, in which the part in the 
constant ci due to the factors is explicitly factored out: 

I A^(a)| < JogDod^N^iBod-^-'hF |/ £ J • • • J , (9.2) 

where (3,p are the parameters in (8.1), (8. 2) with Do, Bo bounded in appendix Al, see 
(A1.2)-^(A1.15), by B, for some constant B (which is not the same as c^M^N^ 4 as in 
(A1.15) only because we have factored out the product of the factors f„ associated with 
each of the h nodes or order 1). 

- h i? 

Note that in general the label v_ in A | „ ( a ) is not equal to the total momentum of 

the tree $ (i.e. A (a) is not a Fourier component, as it depends explicitly on a). 
It is not difficult to see (by working out a special example, see [GGM3] for a trivial 
example; furthermore the proof of the theorem in §10 of [CG], although incorrect as a 
proof of the theorem, nevertheless provides a much less trivial example) that the above 
formula and bounds are optimal, so that a bound like: 

| A|J < B h rf li 'e~^-'-\ (9.3) 
for the i^-th Fourier component of A|(a) would be incorrect. 6 " 

Looking at (9.1) we can therefore bound the sum over v_ and h of the second addend 
by first truncating the series to order h < N: the remainder is estimated, as in §8, via 
the first of (8.1), 66 and will be of the order 0(('q- Q ^i) N ) so that if < r)® +1 it has size 
not exceeding 0(rf ), much smaller than the generically expected contribution from 
the first order terms. 

The first N orders can be bounded (as in (8.6) and using d = 772) by: 

JogDo V - f3/2 N i f2^ h F h e-^^i^BoN^-^^hF < 

h=2 iL 1 +-+K h = iL (9-4) 

< 2JogDoV- f3/2 c l N i ^FB N l +P+H\r 1 -^ 2 ) 2 e- K ^\ , 

where we have used that T, „ 1+ ...+ „ h = K ^ < cV^W* for some c, because 

the sum has dimension < hi (in fact = h(l — 1)) and we have assumed that /1 is so small 
that (fiFBoN^P+Hl^-^/ 2 ) 2 < |. Hence under the condition: 

H/V 1 = f | M | 2D o c l N 4 J o gv- 0/2 [fBoN^+H^-^ 2 ^ ' < 1 (9.5) 

The above bound, after summation over h, is claimed to hold in [RW], theorem 2.1, for a quantity that 
the authors appeared to think, in private correspondence, to be quite simply related to the A ^ (a) of 

[G3] and of this note. 

Which descends, for instance, from the KAM— type theorems in §5 of [CG], hence it holds for analytic 
perturbations and it is uniform in the perturbation as long as that is analytic in llma^l < k and 

k is a prefixed quantity, no matter how small. 
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the sum in (9.4) is bounded by: /j,e K l - 1 . Since A j v is essentially 6c e K ' - ' , we can write 

A T (a)= M Y, e^e-*l^l(A^ + fie^T „(aj) 

0.1 t O;\v\<N (9.6) 

+ 0(^ 1 ) + 0( M 2 e -^ 1 ), |U«)|<1, 

where 0(rp 1 ) is a bound on the terms with h > N and 0(fi 2 e^ ! i' n 1 ) is a bound on the 
terms with h < N with a bound proportional to e - ^'-'-' w ith \v\ > N. 
By using the expression (9.6) for the splitting vector we can compute the determinant 
of the splitting matrix or, similarly to [DGJS], the variation of the pendulum energy 
evaluated on the stable and unstable manifolds, both at ip = n. 

In both cases, by conforming to the traditional terminology, we call Mel'nikov integral 
the first (non trivial) order contribution. 

If one computes the determinant of the splitting matrix, one has to differentiate once 
with repect to a the splitting vector A|(a). Since the derivative of the terms arising 
from the first order contributions is trivial we only need to say that the derivatives of 
£ „ ( a ) can be computed graph by graph and in each graph they amount to a multi- 
plication by a component of the total momentum. Therefore under a slightly stronger 
requirement on [i (since \u\ < n^ 1 one has to require that fi is smaller than n times the 
Ho defined above: \fi\ < ^n) one finds 

det9 Q A T (a) = -2(^Y V ^ , • * - r 

- _T ^ = ° { 9 J ^^^sinhffo;.^-!] sinh [$ w • v'g~^ 

■ (jvfv>{v A v'f + Me- K ^l+^'l)r' £;£ ,) + 0(p 2 e-% T 1 ) , 

(9.7) 

with ir'j, v ,\ < 1, under the stated condition on \fj,\ (i.e. \fj,\ < v/io). 
As to the variation of the pendulum energy, by calling M(a) the Mel'nikov integral 
(depending on a, contrary to the homoclinic determinant (9.7), which is computed at 
the homoclinic point ip — ir, a = 0), one has 

M(a)=fi Y e'"-^ 1 -- 1 m ■ &}„ (9-8) 

vjiO;\v\<N 

so that, for some constant C, 

\HL ■ A T (a)-M(a)| < C|m| 2 e"* 1 ^ 1 1 u | e - K ^l (9.9) 

v ^0;| v\<N 

which can be compared with the results in [DGJS], and, unlike [DGJS], it holds also in 
the anisochronous case. 

Note that so far no use has been made of the assumption that I = 3 nor that lj has a 
rotation number equal to the golden mean nor of the contents of the paper [DGJS] . The 
above (9.8) and (9.7) also seem well suited for extensions to I > 3 (or to more general 
rotation numbers) of the Mel'nikov integral dominance. 

Remark. Assuming f E = e~ K l-l for \v\ < N, by an argument of [DGJS] (exposed in 
§6 of [DGJS], see also §2) we see that, if I = 3 and co has rotation number equal to 
the golden mean (wqi/wo2 = |(\/5 — 1)), in (9.8) the Mel'nikov integral dominates the 



6c Precisely A* = -iu} v 2ZL . -j^-= r = -iu f V* v e"f» 'l^l with 1 < < (2| w ■ 

— \v_ — — y sinh -Sg _1 u> ■ v\ — — — — — 



a\g- L + i). 
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asymptotics as r\ — > if /i is small as indicated. In fact, apart for some "special" small 
intervals of values of r\ and of a, only two addends matter in (9.8) (see [DGJS]): they 
correspond to v equal to ±v_ k — ±(fk+i,—fk) where fj is Fibonacci's sequence and 
k = k{r}) is the value that minimizes 

Kkj-l + fT^Wcr (9-10) 

(of course \vk\ < N): the corresponding miminum has the form c(n)n~i, with c(rj) > an 
explicitly computable quantity bounded above and below by positive constants, uniformly 
in r\ — > 0. The special values of n are those for which the minimum is taken over two 
values of fc, or of 77's close enough to such values; the special values of a are those for 
which | sin a-v k \ < b, where b > is an arbitarily prefixed quantity (if b is too small then 
the size of n for which dominance occurs becomes correspondingly smaller). Excluding 
the just mentioned values of a and r\ the asymptotics of the Mel'nikov integral M(a) 

is proportional to /j,rj~i e~ c ^ v 4 sin v_ k ■ a. Moreover to the r.h.s. of (9.9) the same 
considerations trivially apply, i.e. only two (or four) addends matter so that it is bounded 

i -- 

by a constant times i^vf 3 e ~ c (v)v 4 and the first order dominance for u- A»(a) follows, 
for all a except the special ones. So [G3] and the first order analysis (in [DGJS] 6d ) prove 
the first order dominance not only in the isochronous done in [DGJS], but also 

in the anisochronous case. 

But of course the really interesting quantity is (9.7): i.e. the actual splitting at a = 0. 
For this quantity we cannot conclude dominance of the first order term, unless some 
(likely, see [GGM1], but yet undiscovered) cancellations take place: it is manifest from 
(9.7) that the would be leading term (v_ = v_' = v_ k ) vanishes. Hence no conlusion can 
be drawn, in the class of models considered here, about the splitting at the homoclinic 
point (besides the one already derived in §8 and saying that it is smaller than any power 
of n). Unless one considers 77 — > along a special sequence in correspondence of which 
there are two pairs of successive Fibonacci's vectors with the same, or almost the same, 
K l^fcl + fl^'^fel : m tms case : obviously, we have dominance of the Mel'nikov term 
because no two Fibonacci's vectors are parallel, see (9.7). Note that the expressions 
"almost the same" and "close enough" , used above for the variable 77, simply mean that 
the two Fibonacci's vectors such that the quantity (3k = k\ v_k\ + f 1^ ' H-k\ i s closest to 
its minimum over k for two values of fc, say fc, fc + 1, and pk — Pk+i is a quantity that is 
uniformly bounded as n — > 0. 

§10. Comments. 

(1) The above theory is different from that of [DGJS] in a radical way; in fact it works in 
the isochronous case (J = +00) as well as in the anisochronous one. This is a new result 
(see however the acknowledgements below) . And the last remark shows that the problem 
of the homoclinic splitting asymptotics with 2 fast forcing frequences is still open if it 
is intended as defined by (9.7) (as one should in view of the possible applications to 
heteroclinic chains and Arnol'd diffusion). This is so in spite of apparent claims to the 
contrary, [St]. 

(2) In three time scales problems (i.e. when lj = (rj~^u>o\, rpwo2) and g = 0(1)) the 
above arguments do not lead to any conclusion because u ■ v_ is not always proportional 
to n~ 2 but it can even be proportional to rj 2 also for small \v\: so that the exponentially 
small factor is not always present. This means that even the individual components of 



The rest of the paper [DGJS] provides an interesting partial alternative to the work in [G3]; the partiality 
being mainly due to the isochrony assumption, essential in the theory of [DGJS] . 
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the splitting AJa) in general are not only not dominated by the Mel'nikov integral but 
have size of the order of a power of rjl 

Regarding the intersection matrix, if (9.3) could be assumed valid three matrix elements 
would be exponentially small and also the determinant would be (exponentially small 
and) dominated by the Mel'nikov integral. But, as we said, (9.3) is false. In fact the 
intersection matrix has, in this case, all matrix elements that have polynomial size as 
77 — > 0, starting with second order: this part of the theorem in §10 of [CG] is, of course, 
valid. 66 And it is perhaps natural to think that also the splitting (9.7) has size of order 
of a power of 77. This is however an error, as the determinant is indeed exponentially 
small in n, although the matrix elements have size of the order of a power of n, due to the 
presence of cancellations in the determinant. This was proved in [GGM1], where such 
cancellations are exploited. 

The above error was in fact suffered in [CG] : a trivial computational error in one of the 
about twenty terms contributing to the third order crept in and led the authors of [CG] 
to believe they had checked that even the determinant of the splitting matrix would be 
polynomially bounded away from zero as r\ — > (i.e. bounded as the individual matrix 
elements are) under the convergence condition that /i is itself a (large) power of 77. In 
any case these considerations show that there is apparently no hope, without exploiting 
cancellations in the determinant, to obtain a exponentially small bound (the matrix 
elements have polynomial size) ; the derivation of the exponentially small bounds for this 
case from (9.3), to which many authors seems apparently to believe, is false because (9.3) 
is false, in general and on the section considered here. 

(3) But the paper [CG] had, nevertheless, laid the foundations of the theory developed in 
[G3]: a very simple and systematic theory, because of its field theoretic viewpoint, and 
very flexible thanks to the possibility of rapid (graphical) comparison of arbitrarily high 
order terms of perturbation expansions. Enough to cope with the error and to correct 
it, as shown in [GGM1] (an achievement scarsely appreciated in [St]). In fact the papers 
[G3], [GGM1] suggest which could be the general solution to the problem of the splitting 
(9.7) even when the interaction is a trigonometric polynomial: see the conjecture in [G4]; 
a "few" cancellations remain to be checked to prove (or to disprove) the conjecture. 

In fact the reason the error was not spotted by the authors of [CG] (it was pointed out 
by Gelfreich) was due to the apparent matching of the result with the intuitive idea that 
the splitting in three time scale problems would be "large" . The splitting, defined as the 
smallest eigenvalue of the intersection matrix, has the dimension of an action (being the 
derivative with respect to an angle of the difference in the actions for the stable and for 
the unstable manifolds) hence it has to be measured with respect to an action. In [CG] 
the relevant action with which to compare the splitting was the size of the gaps, in action 
space, between the average actions of the invariant tori surviving the perturbation. 

The latter gap was estimated in [CG] by a power, as large as wished, of 77 times 
max{Jog, Jg}'- hence the splitting was larger than any power of v^ 1 with respect to 
the gaps. 

The proof was wrong and it was corrected in [GGM1]. And in [GGM2] it has been shown 
(to the skeptics, see [St]) 6 ^ that the gaps have in this case size of almost 0(e°^~ ri >) 

while the splitting is dominated by the Mel'nikov integral for generic perturbations, even 

_ 1 ^ w 

if of polynomial type: the Mel'nikov integral being of order 0(e~ v ~ ) if the frequency 
vector is lj = (n^ui,n^^uj). So the splitting is still very large. 69 



The theorem was "proved" by showing that three matrix elements had polynomial size while the fourth 
was "exponentially small" : the last claim is incorrect because of a trivial computational error which, if 
corrected, yields that even the fourth matrix elements is of polynomial size. 

( 't Extending a property and a method that was noted long ago by Neishtadt, [Ne]. 

b 9 And even larger than value claimed in [CG]: the ratio between the splitting and the gap being now 
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(4) Recently we learnt of Eliasson's remark that the homoclinic splitting Q ( a ) is a 
gradient Q (a) = 9 a $(a) of some potential $: this is due to the Lagrangian nature 
of the stable and unstable manifolds. Once this is known one recognizes that the tree 
expansion in sections l-j-8 of [G3] provides a power series expansion of the potential 
see [G4]. 

(5) One can apply the above theory to the case in which the dimension of the angles a 
is 1: this is the periodic forcing case discussed in §8. In this case there is no distinction 
between the (9.7) and (9.8) and one finds an asymptotic expansion for such splittings (see 
§8). This case is closely related to a corresponding problem in the theory of the standard 
map: the basic approach is due to Lazutkin and recently there has been renewed interest 
in it. The paper [Ge], which is nice and accurate, also fills some (apparently) missing 
points in the original work. The approach in [Ge] is also based on formal power series 
expansions for the separatrix; but unlike [G3] the control of the remainders is not solely 
made via the analysis of formal series, thus providing an interesting alternative, although 
(of course) the amount of work is comparable. 

The drawback, at the moment, is the need of an "integrable normal form" near the 
invariant tori (circles in this case). Therefore, although adapting the method of [Ge] to 
cover the cases of Hamiltonian systems appears straightforward in the case of isochronous 
systems (which trivially admit an integrable normal form), a possible extension to the 
anisochronous case seems to require substantial extra work, at least if one proceeds along 
the lines of [Ge]. The method of [G3] has not (yet) been applied to study the standard 
map case. 

Likewise the extension to higher dimension seems to require further extra work (a begin- 
ning of the theory is in [DGJS]); and so does the extension to three time scales problems. 

Once the above extensions will exist they will be be very interesting, and, solving in a 
substantially different way the problems solved in the latter papers, they will permit a 
(better) comparison with the methods in [G3] and [GGM1]. 

(6) Finally the reason why one may hope to get large absolute 611 splitting in three degrees 
of freedom systems is not because some uj ■ v_ (in the small factors e - ' — '-'2?) becomes 
small when v_ becomes a good rational approximation of uj as sometimes hinted in the 
literature, see [St]. This is shown by the three time scales problem in (2) above in which 
lo ■ v_ can be very small even for small v_ (i.e. if the fast component of y_ vanishes) and 
nevertheless the splitting does not have a polynomial absolute size, due to remarkable 
cancellations (see [GGM1]). One has rather to think that the splitting can become large 
because the perturbation f(a,ip) has a small analyticity band in complex ip plane: in 
fact the higher the degree No in cp the smaller it seems that one has to choose /x with 
respect to rj if one wants to guarantee Mel'nikov's approximation to dominate. 
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